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ABSTRACT 

Microwave  attenuation  in  rain  at  C-band.  ob¬ 
served  in  radar  data  from  Alberta,  Canada,  has 
been  assessed  by  means  of  two  different  estim¬ 
ators:  one  employing  Kdp  measured  simultan¬ 
eously  by  an  S-band  radar  and  the  other  a 
comparison  of  synchronised  S-  and  C-band  re¬ 
flectivity  measurements.  While  good  agreement 
between  the  two  has  been  found  in  some  cases, 
in  others  the  second  method  gave  a  larger  estim¬ 
ate  than  the  first  by  up  to  a  factor  of  two.  It 
is  possible  that  the  Kdp  method  is  underestim¬ 
ating  the  attenuation  due  to  the  assumption  of 
a  Marshall-Palmer  drop  size  distribution  which 
does  not  accurately  enough  represent  the  true 
drop  size  distribution  present. 
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1  INTRODUCTION 

Knowledge  of  microwave  attenuation  in  rain  is 
important  not  only  for  predicting  communica¬ 
tion  link  performance,  but  also  for  interpreting 
radar  images.  In  particular,  the  role  of  atten¬ 
uation  at  C-band  is  important  for  interpreting 
data  from  many  operational  radars.  This  pa¬ 
per  aims  to  investigate  C-band  attenuation  by 
analysis  of  data  from  two  co-located,  synchron¬ 
ized  radars  in  Alberta,  Canada.  One  of  these 
radars  is  a  conventional,  linearly  polarised  C- 
band  (5.635  GHz)  radar  while  the  other  is  an 
S-band  (2.88  GHz)  radar  employing  circular  po¬ 
larization  diversity.  If  the  C-band  attenuation 
can  be  successfully  estimated  from  the  polarisa¬ 
tion  information  at  S-band,  in  particular  Kdp^ 
then  it  should  be  possible  to  extend  the  approach 
to  estimate  attenuation  at  other  frequencies. 


2  ESTIMATION  OF  C-BAND  ATTEN¬ 
UATION 

In  a  region  filled  with  rain,  with  a  drop  size  dis¬ 
tribution  iV(D),  horizontally  (H)  and  vertically 
(U)  polarized  waves  have  propagation  constants 
given  by 

KH.v  =  +  fHy{D)N{D)dD  (1) 

where  ibo  =  X’  ^  space 

wavelength  of  the  radar,  and  fjjy  arc  the  ap¬ 
propriate  forward  scattering  amplitudes.  For 
ground-based  radars  and  raindrops  of  oblate 
spheroidal  shape,  as  is  thought  to  be  the  case 
for  all  but  the  smallest  drops,  Kfj  and  kv  will  be 
different.  This  leads  to  both  a  differential  atten¬ 
uation  and  a  differential  propagation  phase  shift. 
It  is  the  latter  which  is  of  interest  here. 

The  two-way  differential  propagation  phase 
shift  0DP  through  a  uniformly  rain-filled  region 
of  length  L  is  given  by 
180 

=  2  X - Re(Kii  -  Kv)  X  L  (2) 

TT 

Unlike  measurements  of  backscatter  quantit¬ 
ies,  such  as  reflectivity,  which  relate  to  a  single 
bin,  0DP  is  accumulated  over  increasing  num¬ 
bers  of  bins.  We  may,  however,  derive  the  quant¬ 
ity  Kdp  which  relates  to  a  single  bin  by  using 

..  ,  ODpU)  -  <^DP(j  -  /o\ 

RdpU)  = - 2A -  ^ 

where  <pDp(j)  refers  to  the  differential  propaga¬ 
tion  phase  shift  along  a  path  from  the  radar  to 
the  jth  bin  and  A  is  the  bin  length.  Kdp  is  the 
specific  one  way  differential  propagation  phase 
shift  and  it  is  from  this  quantity  that  the  C- 
band  attenuation  has  been  estimated  using  an 
empirical  relationship  derived  by  Tan  [1].  This 
relationship  is 

UH  =  0.248A'^|>^'  (dB/km)  (4) 

where  a/f  is  the  specific  C-band  horizontal  at¬ 
tenuation  and  Kdp  is  at  S-band.  A  Marshall- 
Palmer  drop  size  distribution  (DSD)  [2],  drop 
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Figure  1:  Percentage  error  in  attenuation  estim¬ 
ate 


shapes  deduced  from  measurements  by  Goddard 
and  Cherry  [3]  and  Pruppacher  and  Pitter  [4] 
and  a  temperature  of  0  C  have  been  assumed 
in  the  derivation  of  this  equation.  The  most 
important  of  these  assumptions  is  that  relating 
to  the  drop  size  distribution  as  this  can  lead  to 
a  significant  error  in  the  attenuation  estimate. 
This  is  illustrated  in  Figure  1  for  which  expo¬ 
nential  drop  size  distributions  of  the  form 

N{D)  -  No  exp  ^  ^  j  (5) 

were  used  to  calculate  directly,  using  the  forward 
scattering  amplitudes,  both  and  Kdp  the 
appropriate  frequencies  for  a  range  of  Do  and 
No  values.  The  calculated  Kdp  was  then  used 
in  equation  (4)  to  find  an  estimate  of  the  attenu¬ 
ation,  aff(estimateh  and  contours  of  the  percent¬ 
age  error  in  the  estimate,  defined  to  be 


(b) 


were  drawn.  Although  the  Marshall  Palmer  DSD 
has  a  value  of  No  equal  to  8000  ,  the 

contour  of  zero  error  is  not  a  vertical  line  through 
this  value  because  the  power  law  of  equation  (4) 
is  not  exact  but  a  least  squares  fit.  It  can  be  seen 
that  an  underestimate  of  the  attenuation  by  50% 
would  not  be  an  unreasonable  occurrence. 

It  should  be  noted  that  in  rain  regions  Kdp 
will  always  be  positive  and  so  this  method  of  es¬ 
timating  attenuation  is  clearly  limited  to  regions 
through  which  the  measured  <i>DP  is  monotonic- 
aily  increasing.  Since  the  measured  Odp  includes 


a  contribution  from  backscatter  phase,  it  can  de¬ 
crease  in  regions  containing,  for  example,  hail  or 
ice  particles. 


3  RESULTS 

To  investigate  the  performance  of  equation  (4)  in 
practice,  a  second  estimate  of  attenuation  was 
obtained  by  comparison  of  the  S-  and  C-band 
reflectivities.  It  was  assumed  that  in  regions  of 
rain,  if  depolarization/attenuation  effects  were 
removed,  the  S-band  circular  reflectivity  should 
not  be  significantly  different  from  that  of  the 
C-band  horizontal  reflectivity.  The  S-band  re¬ 
flectivity  was  corrected  for  propagation  effects 
using  the  method  described  in  [5],  Where  C- 
band  reflectivity  was  observed  to  decrease  relat¬ 
ive  to  the  S-band  reflectivity  over  a  number  of 
bins,  this  decrease  was  attributed  to  the  effects 
of  attenuation.  Regions  of  rain  and,  in  partic¬ 
ular,  those  where  severe  attenuation  might  be 
observed  were  found  by  identifying  regions  of 
rapidly  increasing  Figure  2  shows  an  ex¬ 

ample  of  such  a  region  from  a  storm  on  21  July 
1979  where  odp  increased  from  9*  to  87  over  10 
bins  or  a  tw*o-way  distance  of  21  km.  This^  gives 
an  average  Kdp  through  this  region  of  3.7  /km, 
which  using  equation  (4)  gives  a  total  two-way 
attenuation  through  this  region  of  16  dB.  A  more 
exact  bin-by-bin  calculation  gives  the  total  at¬ 
tenuation  to  be  15.1  dB,  which  compares  very 
well  with  the  decrease  seen  in  C-band  reflectiv¬ 
ity  relative  to  that  at  S-band  of  14.9  dB.  Figure 
2  also  shows  the  S-band  reflectivity  after  correc¬ 
tion  for  propagation  effects,  the  original  C-band 
reflectivity  and  the  C-band  reflectivity  corrected 
bin- by-bin  using  equation  (4),  Apart  from  two 
bins  where  the  C-band  reflectivity  has  been  over¬ 
corrected  by  about  6  dBZ  the  two  are  in  excellent 
agreement. 

Such  good  agreement  is  not  always  seen  as 
is  illustrated  in  Figure  3  which  uses  data  from 
storms  occurring  on  6  different  days  during  the 
summer  of  1979.  Each  data  point  shows  the  at¬ 
tenuation  estimated  using  equation  (4)  bin-by- 
bin  through  a  region  of  increasing  <Pdp  plotted 
against  that  found  by  comparison  of  reflectivit¬ 
ies.  The  solid  line  indicates  where  the  two  estim¬ 
ates  are  in  agreement.  The  minimum  number  of 
bins  over  which  the  attenuation  has  accumulated 
is  5,  while  the  maximum  is  14.  The  maximum 
attenuation  estimated  by  equation  (4)  was  18.9 
dB  which  took  place  over  12  bins,  while  the  max¬ 
imum  estimated  from  differences  in  reflectivity 
was  39.6  dB  which  took  place  over  11  bins. 

In  many  cases  equation  (4)  underestimates  the 
attenuation  with  respect  to  that  found  from  dif- 


100.0 


80.0 


60.0 


40.0  !- 


20.0 


C-band  reflectivity  corrected  (dBZ) 
S-band  reflecdvrty  (dBZ) 

C-band  reflectivity  (dBZ)  ,  ' 

Differential  phase  shift  (degc^^^  \ 

/  \ 


0.0 


30.0  35.0  40,0  45.0 

Bin  number 


50.0 


Figure  2:  Data  from  storm  of  21  July  1979 


Figure  3:  Comparison  of  two  estimators  of  at¬ 
tenuation 


ferences  in  the  reflectivities.  In  view  of  Figure  1, 
however,  this  could  be  attributed  to  the  fact  that 
the  DSDs  are  not,  in  reality,  those  of  Marshall- 
Palmer. 

Although  in  1979  the  polarization  parameter 
data  was  not  considered  reliable  enough  to  use 
in  estimating  drop  size  distributions,  data  from 
1991  was  available  which  could  be  used  for  this 
purpose.  A  DSD  of  the  form  given  in  equa¬ 
tion  (5),  where  A^o  and  Do  were  free  to  take  any 
values,  was  assumed.  The  parameter 

Zdr=  (7) 

which  can  be  derived  from  the  circularly  polar¬ 
ised  data  is  independent  of  No  and  was  used  to 
determine  the  value  of  Dq.  Nq  was  then  de¬ 
termined  using  either  Z  or  A'/jp.  It  was  found, 
however,  that  in  many  cases  the  estimate  of 
No  from  Z  did  not  agree  with  that  from  Adp* 
This  lack  of  agreement  is  illustrated  in  Figure 
4  where  ZzfB  is  plotted  against  Ze|Kop^  where 
Zfi,  the  effective  reflectivity  factor,  is  in  units  of 
Both  these  quantities  are  independ¬ 
ent  of  iVo  and  a  theoretical  curve  calculated  as¬ 
suming  an  exponential  DSD  and  a  range  of  Do 
from  0.5  to  5  mm  has  been  plotted  together  with 
measured  values  from  a  storm  on  26  July  1991. 
The  measured  values  are  all  from  regions  of  the 
storm  through  which  (j)DP  was  increasing.  Also 
plotted  are  theoretical  curves  for  gamma  DSDs 
of  the  form 

N(D)  =  D^No  exp  ~  ) 

for  /x=2  and  ^=5.  It  can  be  seen  that  the  curve 
is  insensitive  to  the  form  of  DSD  assumed  and 
yet  many  of  the  points  lie  very  far  from  it,  with 
the  values  of  Z^l K dp  being  much  higher  than 
expected.  One  possible  reason  for  this  is  that 
the  region  may  not  be  entirely  composed  of  rain 
but  may  include  hail  or  melting  hail.  Further 
investigation  of  this  is  required.  The  presence 
of  hail  could  also  explain  the  underestimate  of 
attenuation  by  Kdp  compared  to  that  given  by 
the  difference  in  reflectivities.  This  is  because, 
after  removing  propagation  effects,  the  S-  and  C- 
band  reflectivities  would  not  necessarily  be  the 
same  in  a  region  containing  hail. 

Figure  5  shows  a  typical  case  where  C-band  re¬ 
flectivity  including  attenuation  effects  has  been 
calculated  using  an  exponential  DSD  derived 
from  (Z,Zdr)  and  from  (Kdp^^dr)-  The  meas¬ 
ured  C-band  reflectivity  is  shown  for  comparison. 
It  can  be  seen  that  using  Z  to  derive  No  rather 
than  Kdp  leads  to  severely  overestimating  the 
attenuation  observed,  while  using  Kdp  improves 
the  situation  but  leads  to  an  underestimate.  Pos¬ 
sibly,  the  presence  of  hail  has  lead  to  the  S-band 


-  Figure  4:  Comparison  of  theory  and  measured 
radar  parameters  for  a  storm  on  26  July  1991 


Figure  5:  Measured  C-band  reflectivity  com¬ 
pared  with  that  predicted  using  drop  size  dis¬ 
tributions  derived  from  S-band  measurements 
(from  storm  of  26  July  1991) 


reflectivity  being  enhanced,  while  Kdp  been 
reduced. 


4  CONCLUSIONS 

Significant  attenuation  has  been  observed  at  C- 
band  by  comparing  the  reflectivities  measured 
by  the  two  co-iocated  radars.  The  empirical  re¬ 
lationship  derived  by  Tan  [1]  appears  to  predict 
attenuation  well  in  some  cases  but  often  may  be 
underestimating  it.  This  may  be  due  to  the  as¬ 
sumption  of  a  drop  size  distribution  of  Marshall 
Palmer  form,  inherent  in  this  equation,  being  an 
inappropriate  choice.  An  attempt  to  verify  this 
hypothesis  was  inconclusive  as  the  three  meas¬ 
ured  quantities  of  Z,  Kdp  and  Zdr  were  not 
consistent  with  a  model  of  either  an  exponential 
or  gamma  drop  size  distributions.  This  may  be 
due  to  the  presence  of  hail  or  melting  hail,  and 
further  investigation  is  required. 
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